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Consider  the  differential  equation  (1)  x  =  f(x)  in  a.  Banach  space  and  let 
x*  be  an  equilibrium.  The  basic  question  treated  is  the  following:  if  x*  is 
asymptotically  stable  and  if  (2)  x^^  =  x^  +  hcp(x^,h)  is  a  one-step  method,  with 
fixed  step  size  h,  for  integrating  (1),  then  does  the  sequence  converge  to 

x*?  It  is  shown  that  uniform  asymptotic  stability  of  (1)  implies  stability  of  (2) 
and  that  exponential  asymptotic  stability  of  (1)  implies  asymptotic  stability  of  (2) 
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Abatraot 

Conaidar  the  differential  aquation  (1)  k  •  f(x)  in  a  Banaoh 
■paoe  Aud  lot  x*  bo  un  equilibrium,  Th«  baaio  question  treated  lo  the 
followin&t  if  x*  i«  uaymptotioally  etable  And  if  (2)  x^^  ■  x^  + 
ttp(x^,h)  if  a  ono-atop  method,  with  fixed  atep  aieo  h,  for  integrating 
(l),  then  dowa  tho  aequenoe  x^  oonvorge  to  x*1  It  i«  «ho\m  that  uniform 
asyaptoMn  at Ability  of  (l)  implies  atAbility  of  (&)  And  thAt  exponential 
aaymptoUc  stability  of  (1)  impliea  asymptotic  stability  of  (2). 
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1.  Introduction, 

Conaid«r  the  differential  equation 

(1.1)  x  »  f(x) 

and  let  x*  be  an  equilibrium  point.  The  basic  question  to  be  treated 
here  is  the  following!  if  x*  is  an  asymptotically  stable  equilibrium  and 
if 

(1.S)  Vi  ■  *k  +  W’<Vh)  " 

is  a  one  step  method,  with  fixed  step  size  h,  for  integrating  (1.1),  then 
does  x^  converge  to  x*  as  k  tends  to  infinity?  We  shall  show  in 
our  first  main  theorem  that  uniform  asymptotic  stability  of  (1.1)  implies 
stability  of  (1.2)  and  in  our  second  main  theorem  that  exponential  asymp¬ 
totic  stability  of  (l.l)  implies  asymptotic  stability  of  (1.2)  (improving  ' 
a  result  of  Skalkina  [ 11]) . 

% 

Our  interest  in  the  problem  considered  here  stemmed  from  an  investi¬ 
gation  of  iterative  methods  for  solving  the  equation  F(x)  =  0  in  a  Banach 
spec*.  If  f(x)  is  a  function  whose  zeros  include  the  zeros  of  F  (for 
exaaple,  f(x)  =  -(F^)  -1F(x)),  then  numerical  integration  of  (l.l)  will 
lead  to  iterates  corresponding  to  points  x 
curve.  If  the  initial  point  x^  is  in  a  region  of  attraction  of  the  equi¬ 
librium  x*,  then  under  what  condition  does  x^  converge  to  x*?  Various 


on  the  solution 


■wt  *--■** 


author*  have  used  similar  idea*  to  develop  algorithms  for  solving  F(x)  »  0 
([  1,2,3, 5#10'J)  in  particular  situations.  For  uxample,  Boggs  ([!])  has  in- 
tegrated  the  equation  jc  »  -(F')  T(x)  with  the  A -stable  methods  of 

A 

Dahlqui at  to  generate  iterates  which  converge  to  a  root  of  F,  In 
[2,3,10],  Euler  and  Runge-Kuttu  integration  methods  are  used  to  generate 
iterates  x^  which  eventually  lie  within  the  region  of  convergence  of 
Newton's  method.  Here,  results  are  devoloped  for  general  one  step  methods. 

2.  Uniform  Asymptotic  Stability. 

I-et  X  be  a  real  Banach  space  with  norm,  ||*||,  and  let  S(r)  «* 

(x|  )|xl|  <  r)  be  the  closed  ball  of  radius  r  about  0  in  X.  We  let  f 
be  a  mapping  of  X  into  itself  and  x*  be  a  zero  of  f.  We  assunie, 

.  .  -  '\ ■  ■  i  .V 

without  loss  of  generality,  that  x*  =  0.  Now,  suppose  that  f  i3  defined 
on  the  ball  S(R)  and  that  <p(x,h)  is  a  mapping  of  S(R)  x  [0,h0l  into 
X.  We  assume  throughout  the  sequel  that  the  following  conditions  are  sat¬ 
isfied! 

(2.1)  there  are  positive  constants  L  and  L’  such  that  ||f(x) -f(y)|j  < 

L|| x-y|!  and  ||cp(x,h) -<p(y,h)J|  <  L’||x-yj|  for  all  x,yGS(R)  and 

°  <  h  <  hQ; 

(2.2)  cp(:t,h)  is  uniformly  continuous  on  S(R)  x  [0,13q]j 

(2.3)  T>(x,0)  =  f(x)  for  all  x  G  S(R);  and 

(2.4)  f(0)  «  <p(0,h)  =  0  for  all  he  [0,hQ]. 
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We  consider  the  differential  equation 

(2.5)  x  ■  f(x) 

and  the  one  step  integration  method 

(2.6)  xk+1  *  xk  +  h?s(xk,h). 

,■  v 

[Note  that  (2,6)  is  consistent  in  viev  of  the  assumption  (2.3)],  We  now 
have 

DEFINITION  2,7.  The  solution  x  -  0  of  (2.5)  is  uniformly  stable  if,  given 
€  >  0,  there  is  a  5(c)  >  0  such  that  ||x^||  <  6(e)  implies  that 
||x(t)t0,XQ)||  <e  for  t  >  tQ  where  xCtjtg,^)  is  the  solution  of  (2.5) 
with  x(tQj  VV  =  X0*  The  solution  x  ■  0  of  (2.5)  is  uniformly  asymp¬ 
totically  stable  on  a  ball  S(r)  if  it  is  uniformly  stable  and  if,  given 

i 

c  >  0,  there  Is  a  T(e)  >  0  such  that  ||x^||  <  r  implies  that 
||x(t>t0,x0)||  <  e  for  t>tQ  +  T(e). 

We  note  that  since  X  may  be  infinite  dimensional,  uniform  asymp¬ 
totic  stability  and  asymptotic  stability  are  not  eq\ii  valent  (( 9])  • 

We  now  assume  that  the  solution  x  =  0  of  (2.5)  is  uniformly 
asymptotically  stable  on  the  ball  S(R)  for  some  R  >  0.  If  6(e)  and 
T(e)  are  the  functions  characterizing  the  stability  of  (2.5)  as  in  defini¬ 
tion  2.7,  then  we  may  assume  that  6(*)  and  T(*)  are  strictly  monotonic 
continuous  functions  (see  [7,  P-  5091).  We  also  suppose  for  simplicity  that 
t0  -  0  and  we  let  x(t>xD)  »  x(t,0,x0).  We  then  have: 


k 


U5MMA  2.8.  let  r,b  be  reel  numbers  such  that  0<r<b<6(R).  22 
there  le  a  tx  >  0  euch  that  inf(||x(tjx0)|||  t  e  (O.tj),  r  <  ll^ll  <  b) 
la  strictly  positive. 

Proof:  Since  b  <  6(R),  II  (|  <  »  and  so  l|x(t;x^)  -*01I“ 

||nf(K(8))-f(x0)}dS  +  tf(x0)||  <  L  /  lloc(s)  -^lids  +  m^)\\  (vhere 
;,o).  It  follows  from  Gronwall' s  inequality  and  an  integration  by 
parts  that 

(2.9)  ||x(t)  -XqII  <I^!^|jte 

Therefore,  ||x(t,^)||  >  llxjl  («*•“)  “d  my  °h°°Se  *1  >  °  SUCh 

^1 

that  1  -  Ltxe  >  0. 

/rni\  „«  i  n(  be  a  continuous  strictly 
Following  Mas  sera  ([9])>  we  -'-et  ' 

,  .  . ..  «/_,\  2r  G(o')  =  0  and  we  introduce  the 

increasing  function  with  G(r)  < 

Lyapunov  function  V(0  for  (2.5)  given  by 

(2.10)  V(^)  -  sep(0(||x(tiXo)||)(l+2t)/a+t)|t>0) 

for  0  <  ||x0l|  <  0  where  p  =  mind,  6(H)). 

LEMft  2.11.  (193).  V(-)  tee  the  following  properties:  (1)  °(II*qID  <■ 
v(jr )  <S,|x0ll;  (11)  llv(*0)-v(y0)ll  f  MH*0-y0"  £2 H-SSSS  M  >  0;  'Ul) 

V(x  )  =  llm  sup  [V(x(k;Xo))-V(x0))/h<  .0(11^11)  (lt2T(ll*oll/2))‘  i  Hi 
0  k  -*  0+ 

(IV)  v(x(k;x0))-v(x0)  <  -M3(!!x0l|){i+2Kia:(!l*(l'-,=0)ll/2))'2  £21 
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ll^liJyoll  <  p  «nd  k  >  0. 

Letting  i|r(||*J|)  =  G(||xJ|){l+2T(||x0||/2)}-2,  we  have* 

LEMMA.  2.12.  If  0  <  r  <  5(p)  and  e  >  0,  then  there  is  a  k(r, e)  >  0 
such  that 


(2.13)  V(x(kjx0))-V(xQ)  <  k{-*(||xJ!)+eG(Hx0||)} 

for  0  <  k  <  k(r, e),  r  <  || xjj  <  8(p). 

Proof;  Choose  t^  >  0  by  lemma  2.8  so  that  m  «  inf(||x(t,*x0)||  |  0  <  t  <  t^, 
r  <  || jo ||  <  8(p)  j  >  0.  Then  0  <  m  <  ||x(t;xQ)||  <  p  for  0  <  t  <  t^  and 

r  <  !!X0H  <  6(p)- 

2 

Since  A(k,  cr)  =  [  l+2k+2T(cr)]  “  is  uniformly  continuous  on  [ 0,  t^]  x 
[m/2,  p/2],  there  is  an  r\  =  rj(e)  such  that  | A(k,o’) -A(0,cr) |  <  e  if 
)k|  +  |  cr  *  -cr|  <  T.  Let  k(r,e)  be  the  smaller  of  t^  and  the  unique  posi¬ 
tive  solution  of  k  +  ^L8(p)keIlk  *  tj.  Letting  o'  =  ||x(kj  x0)||/2  and 

0  "  !lx0il/2»  it  follows  from  (2.9)  that  |  u'  -o|  <  ||x(k;^)  -x^U  /2  < 

Lk 

(L6(p)ke  ) /2  and  hence,  by  virtue  of  lemma  2.11,  that  V(x(k; Xq))  -V(Xq)  < 
-lsO(l|x^U)A(le, a*)  <  kC-^(Hx^JI)  +  eG(||xg||)} . 

LEMMA  2.l4.  There  is  an  U  >  0  such  that,  if  x(t>x^)  and  xQ+hf(3^)  are 
elements  of  S(R)  for  0  <  t  <  h  <  then  H  a^+hfCx^)  -x(hja<^>)||  < 

m>X»- 

Proof:  Apply  Gronwall’s  inequality. 

LEMMA  2.15.  Let  p  »  min{l,  B(R)}  and  suppose  t>iat  <p(x,h)  is  uniformly 
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continuous  on  S(p)  X  [0,11q1-  ££  0  <  r  <  &(p),  th — 

such  that 


there  is  an  h-^r)  >  0 


(2.16) 


Vfjtg+ttPOtyh))  -  V(Xq)  <  -  5  <  0 


whenever  0  <  h  <  h^r),  r  <  U*q11  5  1 


proof,  Assume  without  loss  of  generality  that  S(p)  <  P.  Then,  If  h  < 
(p-8(p))/[6(p)max(Xi,L,)l,  "  el“CntB  °f 

S(p).  Now,  V(VW(V»))-V(V  ilvtV^V-^  -*V**b»l  + 

|V(Vhf(x0))-V(*(hi,0))l  ♦  v(m(h,V)  -  V(^-  «  the  pre- 

,  4V,n4.  -Fnr  0  <  r  <  8(p)  “ul  e  >  0,  there  is  a  k(r,e)  >  0 

vious  lemmas,  that,  for  u  r  \ 

such  that  if 


(2.17) 


<  h  <  h»  .  mln(h(r,e)',h0,(p-6(p))/Kp)'®'tI'.L'l^ 


then 

(2.18)  V(V»(Vt))-,(V  5  Mhll<P(Vh,-f(x0)l1  +  ”n’2!|X;o11 

-  h’Kl!x0ll)  +  «°(lla^Dh 


for  r  <  IIxqII  <  5(p). 

'let  o(h).  »uv(Wx,h)Jp(*,0)»|r<W<M)  and  tah.  «  < 

rtWoWp)*.  ^*’0)  *  ftxM  Sln,‘  "(X'b)  1,u“lfOTly 
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continuous,  a(h)  is  continuous.  Moreover,  a(o)  =»  0.  Thus,  the  equation 
a(h)  +  Nh&(p)  =  i(r)/(4M)  has  a  least  positive  root  h  >  0.  If  0  <  h,(r)  < 
inin(h*,h),  then  it  follows  that 


(2.19) 


V(x0+hcp(x0,h))  -  vC^)  <  -  |ht(r) 


for  0  <  h  <  h^(  r)  and  r  <  ]|  xQ||  <  b(  p) . 

We  can  now  prove  the  following: 

THEOREM  2.20.  Suppose  that  the  solution  x  =  0  of  (2.5)  is  uniformly 
asymptotically  stable  on  S(R).  Then,  for  any  €  >  0,  there  are  h(e)  >  0 
and  K(e)  >  0  such  that  if  ||x^||  <  G(6(p))/2  and  0  <  h  <  h(c),  then  the 
solution  of  (2,6)  starting  from  Xp  satisfies  the  inequalities  (i) 
llxjl  <  p  for  all  k  >  0;  and  (ii)  IjxJI  <  e  for  all  k  >  K(e) /h. 

Proof:  We  may  assume  that  0  <  €  <  6(p).  Let  r  =  G(e)  fk  and  let  h(e)  = 
tnin(h1(r),l/L')  where  h^r)  is  given  by  lemma  2.35.  Also,  let  K(e)  =  . 

2(G(B(p))-G(r/2))/i|r(r). 

We  consider  three  cases,  namely:  (i)  0  <  |) 3C^|}  <  r,  (ii)  r  <  ||x^|]  <  2r, 
and,  (iii)  2r  <  ||x0||  <  G(6(p))/2. 

Case  (i):  If  ||XjJ|  <  r  for  all  k  >  0,  then  (IxJI  <  G(e)  /k  <  e/2  <  e  for 
all  k  >  0.  On  the  other  hand,  if  jlxj!  <  r  for  k  -  0,1, ...,n-l  and 
llxJI  >  r,  then  IjxJ  =  !lxn.p+llcP(xn_1>h)ii  <  HxnJ(1+hL')  <  2r  and  we  re¬ 
gard  xn  as  an  initial  point  for  case  (ii) . 


I 


ft 


Case  ( 1 1 ) i  We  claim  that  ||x^i|  v  «  for  all  k  >  0,  (Note*  that  «  * 

6(p)  <  p,  ]  Clearly  ||>^||  •»*  2r  <  «,  If  r  <  ||x^(|  <«  for  0  k  <  ft, 

then  0(||xntl!|)  <  V(xn+l)  -  V(x0)  ♦  t  Mx^j)  •  Vfx^l  <  V(j^)  . 
(n+l)W(v)/c  *.  v(^)  <  ;?||s^!|  <  l*r  ~  Q(«)  by  virtu#  of  Itmmaa  3,  11  and  6,15, 
Since  G  in  strictly  monotone,  Hx^H  <  *  and  tha  olaim  it  tttabiithtd 
by  Induction, 

Thus,  combining  cast#  (i)  and  (U),  wo  have  show  that  if  ||Xg||  <  0r, 
then  ||xj|  <  f  for  all  k  >  0, 

Case  (iii):  Clearly  HxJI  <  0 (0(p))/2  <  6(p)  <  p.  Suppose  that  r  ;£ 
ll*JI  <  B(p)  for  k  <  n,  Then,  Oflls^J)  <  V(xr^1)  *  V(x^)  +  Wx^j). 

v(xk))  <  vCxq)  -  (n+l)h^(r)/2  <  2)|j^||  <  0(5(p))  by  virtue  of  lamias 
2,11  and  2.15.  Since  G  is  strictly  monotone,  ||x.+ J  <  6(p)  <  p  and  to, 
||xj  <  p  for  all.  k  >  0,  Furthemoro,  if  (n+l)h  >  K(«),  then 
°^*n+Jp  5  v(*0>  -  K(Oi(r)/&  <  £|| XqII  -  Kf«)*(r)/8  <  Q(6(p))  -  K(«)t(r)/»  • 
G(r/2).  it  follows  that  Ijx^JI  <  r/2.  The  theorem  then  follows  frets  tho  t 
first  two  cases. 

We  note  that  the  theorem  does  not  assort  that  the  solution  ■  0 
of  (2.6)  is  ti table  for  fixed  h.  In  other  words,  we  do  not  claim  that 

i 

for  given  h  and  any  e  >  0  there  is  wn  t)«  rj(c,h)  such  that  if 
ll^ll  <  n,  then  ||x^||  <  e  for  all  k.  Bearing  this  in  mind,  we  oonsider 
the  following  two-dimensional  system: 

x  «  y  -  x(x  +  y  ) 


(2.21) 


y  »■*  -  y(xr'  t-  y2) . 


!•#%  V(H|y)  »  »** ,  Then  $(*(t),y(t))  *  .8(x^(t)*ye(t))ft  alarm  80lu. 
tlone  of  (1*11)  end  80)  the  tv i viol  aoluhlen  ia  unifomXv  aaymptotloeliy 
•table,  If  Ruler* a  Mthod  U  applied  to  (B.fcl),  than  tho  difference  eya 


tiia 

(•alt) 


Vi  *  N  *  *»  •  VW 
Vl  ’  •  *“V  ’  VW 


!•  obtained.  Let  ^V(X)  y)  bo  given  by 

<».a»  yfc*)  .  hWu«V  . 


•o  tbot  V(xRtl,y^1)  .  V(j^,yn)  »  ^V(xn,yu).  Ualng  (8.83),  it  la  *uy 
to  varify  that  tha  trivial  aolution  of  (8.88)  la  not  atabla  and  that  all  eolu. 
tluna  with  0  <  h(^yj)  <  1  ♦  (l-h8)1^  (0  <  h  <  1)  ara  attraotad  to  tha 
Invariant  aat  x8+y8  ■  (i-d-h8)1^8)  /h.  Althoufh  tha  trivial  aolution  of  . 
(8,88)  la  not  atabla  for  fixed  h,  tha  aolutlona  of  (8,88)  can  ba  made  to 
remain  arbitrarily  oloaa  to  aaro  by  Initially  ohooalng  h  email  enough. 

In  other  worde,  tha  theorem  aaaarta  that  for  given  «  >  0,  there  la  an  h(«) 
•uoh  that  if  h  <  h(«),  than  the  aolutlona  of  (8.88)  will  11a  within  tha 
ball  8(e), 


3,  exponential  Aaymptotio  Stability. 

We  now  oonaider  the  oaae  of  exponential  aaymptotio  etabillty. 

DEFINITION  3,1,  The  aolution  x  «  0  of  (2.3)  ia  exponentially  asymptotically 
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stable  cm  8(r)  if  there  are  positive  oonntanta  «,  M  suoh  that 

£21  ll^ll  5  *  m»2  !&!*»*&>  th> 

wolution  x^  m  0  of  (0,6)  if  axnonontially  asymptotically  stable  if  there 
are  positive  constants  b,  h^,  Mj,  R  iuoh  that  ||x^||  <  NjJlxJIt"^  for 
oil  k  *  0  whsnavcr  0  <  h  <  »nJ  ||j^||  <  b^, 

Skalkina  (( 11))  ho*  shown  that  if  the  taro  solution  of  (2.5)  1* 
•xponsntlaliy  aaymptotloally  atablt,  than  ao  ia  tha  taro  solution  of  (8.6), 
We  shall  shortly  present  an  improved  version  of  his  result. 

LEMMA  5.8,  If  tha  solution  x  ■  0  of  (fi,5)  is  exponent ially  asynptoti- 
oally  stable  on  S(R),  than  tha  function  W(*)  defined  by 

(3.5)  W(x0)  ■  sup(||x(t|x0)jiexp(arotan  at)Jt>0) 

for  Xq€  8(R)  has  tha  following  properties!  (\)  Hx^H  <  W(xq)  <M||xq||; 

(ii)  | W(x^) -W(yQ) I  <^l*Q-y0lb  Uii)  <  -o’W^)*  and,  (iv) 

i 

wMhjXg))  -  W(xq)  <  -o||Xq|!>\  for  suitable  positive  oonstants  KjO^o' 

(where  a,  M  are  the  oonstants  involved  in  che  definition  of  exponential 
asymptotic  stability) . 

t 

Prooft  Argue  as  in  [7j  pp.  309-511]. 

We  now  have 

THEORKM  3.^.  Suppose  that  the  solution  x  j  0  of  (2.5)  is  exponentially 
asymptotically  stable  on  S(R) .  Assume  also  that  either  (a)  ?(x,h)  « 
f(x)  or  (b)  tp(x,h)  is  (Prochct)  differentiable  in  x  and  cpx(x,h)  i_s 

Uv.r  lol  uniformly  continuous  on  Ffa)  xfO/n^].  Then  the  solution 
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\  ■  o  of  (2.6)  1«  exponentially  asynuatotioally  stable. 

Pvoofi  tot  b«(0,R)  and  lot  K  -  min(Lmllog  R/b,  ((R/b)-l)/laM(L#f),^). 

If  0  <  h  <  K,  then  x  <•  hp(x,h)  and  x  +  hf(x)  are  in  8(R)  for  ||xj|  <  b 
and  ||  x(t,  Xq)  H  <  R  for  all  t>0  if  ||^||  <  b  (a*  ||x(t,j^)||  < 

Kilt1*). 

Now,  lat  a(h)  r  o  or  oup(||9x(x,h)  -  9x(x,C)|||||x||  <  b)  aooording 
ai  hypothesis  (a)  or  (b)  holds.  If  hypothesis  (a)  holds,  then  |  W(x+hp(x,h))  - 
W(x+hf(x))|  ■  0  <  Kha(h)||xj|,  On  the  other  hand,  if  hypothesis  (b)  holds, 
then  Jw(xfhcp(x,h))  -  W(xfhf(x))|  <. Kh||«p(x,h)  -  <p(x,0)||  <  Khj|cp(x,h)  .  9(0, h)  - 
<p(x,0)  +  <p(0,0)||  <  Kh  /  (||9x(tx,h)  -  cpx^tx,0)!!||x||)dt  <  Kha(h)||x||,  [Note 

that  9(0,h)  ■  9 (0,0)  ■  0,]  In  other  words,  we  always  have 
(3.5)  | W(xfhp(x,h))  -  W(x*hf(x))|  <  Kha(h)||xj| 

for  ||x||  <  b. 

Let  h'  be  the  least  positive  root  uf  K[a(h)  +  Nh]  »  c/2  and  let 
hx  be  any  positive  number  with  <  min(h,b‘,2M/c) .  If  h  <  ^  and 
ll^ll  <b,  then  W(^+hp(xQ,h))  -  W(x^)  <  |  WCxQ+hp'x^h))  -  N^hfC^))!  + 

|  W(x^+hf (xQ) )  -  W(x(h;Xg))|  +  W(x(h;x0))  -  W(x^)  <  hJlxJ!  {Ka(h)+KNh-c}  < 

-4 *qII  c/2. 

Now  let  »  M  and  0  *  c/(2l$.  We  ohall  show  by  induction  that 
if  HxJI  <  b M,  then 

-Pkh 


(3.6) 


Hxjjl  <Mj|x^|!e 


for  all  k.  Clearly  (5.6)  hold,  for  k  -  0  and  so,  we  suppose  it  holds 
i  for  0  <  k  <  n.  For  any  euoh  k,  llxjl  <  <  **Jl  <  b  ^  a0» 

1  W(x^1)  -  W(xk)  <  -hllxj.il o/2.  Sinoe  W(xR)  <^1^11,  w(*k+l)  - 

If  W(xK)  -  0  for  any  k  <  n,  then  V(\+t)  •  0  for  &11  and 

i  (3.6)  is  satisfied.  Otherwise,  Ha^JI  <W(xn+i)  <W(xn)(l-Ph)  < 

!  '  w(x0)(l-ph)n+1  <  Mil ^qII  e’ph^n+1\  Thus,  the  theorem  is  established. 

i 

i 


i 
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